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In this letter, the 6- vertex model on dynamical random lattices is defined via a matrix 
model and rewritten (following I. Rostov) as a deformation of the 0(2) model. In the 
large planar limit, an exact solution is found at criticality. The critical exponents of 
the model are determined; they vary continously along the critical line. The vicinity of 
the latter is explored, which confirms that we have a line of c = 1 conformal field theories 
coupled to gravity. 
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1. Introduction 



Starting from tlie observation that the combinatorial properties of large matrix 
models allow them to reproduce summations over discretized surfaces it is possible 
to study various 2D statistical models on random dynamical lattices, that is consider 
systems in which both the "spin" degrees of freedom sitting on the lattice and the lattice 
itself are allowed to fluctuate. Among the models that were solved this way, let us cite 
the Ising model @], the 0{n) model the Potts model AH these models have 

critical points which correspond to conformal field theories of central charge c less or equal 
to 1 coupled to gravity, the limiting case c = 1 being of particular interest. One well- 
known statistical model has until now resisted attempts at an exact solution (on random 
lattices): the 8- vertex model and its critical version, the 6- vertex model. On a flat lattice, 
the 6-vertex model displays an infra-red behavior which spans a whole semi-infinite line of 



c = 1 theories; a similar behavior is expected when put on random lattices IT^. In fact 



a special 2-parameter slice of the 8- vertex model was solved exactly in |T^; the 6-vertex 
point of the model was shown to indeed exhibit a c = 1 behavior. 

In this letter, we shall perform a study of a matrix model which describes the 6-vertex 
model. We shall be concerned with its large limit, which correspond to selecting the 
spherical topology for the lattices. We shall solve the model exactly when the renormalized 
cosmological constant vanishes, that is in the limit where the average size of the graphs 
goes to infinity, and give explicit expressions for some averages. We shall then show how 
to explore the vicinity of this critical region and compute the first correction, which yields 
the string suceptibility. 



2. The model and its saddle point equation 

The model is defined by the following partition function: 

Z = j dXdX^ exp A^tr (^-XX^ + bX'^X^^ + ^{XX^f^ (2.1) 

where X is a general N x N complex matrix. The Feynman rules of the model reproduce 
the configurations of the six-vertex model (figure 1) on a random four-valent lattice. 
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Fig. 1: Part of a diagram generated by the perturbative expansion of (2.1) 



In the parametrization of e.g. |]T^, c/h = 2 cos(A7r/2); the parameter A can, on random 
surfaces, vary in the range < A < 2. For A = we recover the usual 0(2) model 0; 
for A = 1 we recover the critical ABAB model Finally, the singular limit A — 2 

corresponds to the three-colouring problem [ p!3| , p^ . When A is fixed, we can still vary 6, 
which plays the role of (bare) cosmo logical constant. For any A we expect that there is a 
value 6crit(A) for which the average size of the graphs diverges. 

In order to solve the model, we use the following trick: we decouple the quartic 
interaction |T^ 



dAdXdX^ exp 



A^tr 



-XX^ - -A^ + VbA{XXU'^''/^ + X^Xe"^^^/^) 



(2.2) 



where A is a hermitean matrix. This new model is a deformation of the 0{2) model. Its 
Feynman rules (figure 2) allow to interpret it as a model of oriented loops in which each 
left/right turn costs uj^^ = e^*^'^/^; this means that each loop, taking into account its two 
possible orientations, contributes a factor of 2 cos(pA7r/4), where p is the number of right 
turns minus the number of left turns performed when going around the loop. For a regular 
infinite lattice, this number is fixed and we recover this way the usual 0{n) model with 
n = 2cos{pXn/4); but this connection between the 6- vertex and the 0{n) models breaks 
down on a random lattice (because of curvature). It only remains at A = 0, of course, 
where we recover, as mentioned previously, the 0{2) model. 
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Fig. 2: Part of a diagram generated by the perturbative expansion of 
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The next step to solve the model is to integrate out X and X't', and shift A by the 
constant 7 = ^ — : 

Z = J dAdet-^ {u;^A + A^u;'^)exp(^-^tr{A-jf^ (2.3) 



In terms of the eigenvalues of ^: 

N 



/n 



In the large A?" limit, the form a saddle point distribution characterized by a con- 
tinuous density pQ{a)da which fills an interval [a, (3] of the real line. It will be convenient 
to make the following change of variable: a = (5e^. Up to an overall constant we find: 



N 



duie 7; TT (2.5) 

1=1 i^j 



One notices that the two-body interaction now only depends on the difference Ui — Uj. 
More explicitly, we are now trying to minimize an action of the form 

J dup{u) V{u) + \JJ dup{u)dvp{v) V2{u - v) (2.6) 
with the density dup{u) = dapo(«), a potential V{u) = ^{Pe^ — 7)^ and an interaction 

V2{u) = -21og(l-e'^) + log(l + a;V) + log(l+a;-V) (2.7) 

The density p{u) satisfies p{u) > and J dup{u) = 1. It is worth remarking the similarity 
of this problem with the determination of the ground state in the presence of magnetic 
field in Bethe Ansatz solvable models. We shall comment on this analogy later. 

On the support of p{u), the minimization of S leads to the saddle point equation: 

K*p(^i) = /3e"(/3e" -7) (2.8) 

where ★ means convolution product and K is the derivative of — V2: 

K{u) = \ (2.9) 

Principal part at tt = is implied. 

We shall now proceed to solve this equation. 
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3. Exact results at criticality 



The analytic problem ( |2.8| ) is well known to be exactly solvable using the Wiener- 
Hopf technique when the support of p{u) is semi-infinite. In our case this corresponds 
to a = 0, which is precisely the critial regime (i.e. when the renormalized cosmological 
constant vanishes and the area of the lattices becomes large). The support of p{u) is then 
[—00, 0] (figure 3). Note that the situation is then analogous to the Bethe Ansatz study of 
the ground state in the presence of magnetic field of 2D integrable models with only one 
chirality, see for example []T^ . 



Fig. 3: Potential V{u). The eigenvalues fill the well; criticality is attained when 
the eigenvalues start "overflowing" at minus infinity. 



We introduce the Fourier transform K{k)\ 

/+00 
-00 



= 471 



sinh(i + |)7rfcsinh(i - J)7rfc 
sinhnk 



(3.1) 



Similarly we have the Fourier transform p{k) of p{u), which can be defined alternatively 
as an average in our model: 

m = ^(ltrAA (3.2) 



13^'' \N 

It is clear that p{k) is an analytic function in the lower half-plane Imk < 0. We now 
Fourier transform the equation (|2.8|) : 

K{k)p{k) = f{k) (3.3) 

where f{k) is a function whose inverse Fourier transform f{u) satisfies f{u) = f3e'^{l3 e'^ — 'y) 
for u < 0. 
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In order to determine p and /, we decompose K{k) as K{k) = K-{k)/ K^{k) where 
K- (resp. K^) is holomorphic in the lower (resp. upper) half-plane. Explicitly, we choose 



(3.4) 



where u± = ^ ± j and e = log + U- log U- . From ( |3.3|) we infer 

K.{k)p(k)=K+{k)f(k) 

The left hand side is a function which is holomorphic in the lower half-plane, and therefore 
the right hand side must also be. It is now an easy exercise to find it starting from 
f{u) = /3(/3e" - 7), w < 0. The result is: 

where Ki and K2 are K+{k) evaluated at k = i and k = 2i: Ki = cos(A7i-/4) ^~'^ 
K2 = sinlx-K /2) ^~'^'^ ■ Dividing by K-{k) gives p{k). Note however that one must impose 
the normalization condition p{k = 0) = 1. This in fact imposes two constraints since 
generically the function p{k) given by (|3.5| ) has a pole at /c = 0. The pole cancellation 
condition reads 

7 = /^^ (3.6) 

while the normalization condition reads 

m = ^P'K2 = 1 (3.7) 

These two conditions determine the critical values of jS and 7. In terms of the original 
coupling constant b, we find: 

_ 1 sin(A7r/4) 1 
""''-32 Xtt/A cos3(A7r/4) ^^'^^ 

Finally we obtain the expression for p{k): 
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The usual correlation functions (-^trA") are simply given, according to (|3.2|), by 
(3'^p{k = —in). However, to extract the universal information from the correlation func- 
tions, it is simpler to consider the singularities of the density. 

The density po{a) has two singularities: one at a = a = and one at a = /3, which 
correspond to m = — oo and u = 0. 

The singularity of po{a) at a = /3 is determined by the behavior of p{k) as /c — oo. 
We find that p{k) /c""^/^, or 

po(a) ~ (a-/3)i/2 (3.10) 

i.e. the usual square root singularity. 

On the other hand, the singularity of po{a) at a = can be inferred from the poles of 
p{k), which lie at 

i Ti 

k = — n = l,2, ... (3.11) 

u± 

For A > 0, the pole closest to the real axis is k = i/u^, which corresponds to a behavior 
p{u) °° qu/u+ g^j^i^ therefore po(fl) ~ a"+ ^ or more explicitly a leading singularity of 
the form 

2-A 

Po(a)~a2+A (3.12) 

Therefore we have found one critical exponent of our model, which depends continuously 
on A. The other poles give subleading terms in the expansion around a = 0. 

In the limiting case A = 0, we find a double pole at k = 2i, which results in a behavior 

p{u) ~ ue'^'^, or 

Po{a) ~ a log a (3.13) 



4. Vicinity of the critical line 

Outside the critical regime, i.e. when the support [ct, (3] of the density p{u) is finite, it 
is not clear how to solve exactly the equation ( |2.8| ). This is completely analogous to the 
Bethe Ansatz equations for 2D integrable models with 2 chiralities (like massive relativistic 
models). However, what one can do is an exact expansion as a — 0. Here we shall follow 



the method of |]T^. We shall only compute the first correction to the calculation of the 
previous section, up to some constants which could be determined by a more careful study. 

Let us denote B = log{(3/a). We are interested in the limit where B is large (in the 
Bethe Ansatz language, the limit where the two chiralities are almost decoupled from each 
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other). Then we know that we can spht our equation ( p.8|) into two equations for the 
regions u ~ —B and u ^ 0. The first equation is after the shift u ^ u + B: 

K ★ p{u) = Vw G [0, B] (4.1) 

We introduce the function g{u) with support [— oo, 0] such that 

K * p{u) = + e-^'''+g{u) Vw G [-oo, B] (4.2) 

We have included the prefactor so that g{u) has a finite hmit when B — > oo, as a 

simple calculation shows (for A = one should replace with Be~^^). The second 

equation is our usual saddle point equation; we can insert in it the first correction: 

Ki<p{u) = Pe^'iPe'' -'y)+e-^/''+g{u + B) G [-oo, 0] (4.3) 

so that the equation is valid for all negative u. Therefore this is again a Wiener-Hopf 
problem. We shall not bother to solve it explicitly; let us simply write down the form of 
the solution at leading order as S oo: 

= " d'whr) - ^7(Fh) ) + ^-""^^-''kk) (4.4) 

When dividing by K-{k) we are again faced with the problem of the behavior at /c = 0, 
which leads to the two conditions: 

7 = + cie-«/-+ + ■ • • 

(4.5) 

1 = -I3'^K2{1 + C2Se-^/"+ + C3e-^/"+ + ■ ■ ■) 

where ci, C2, C3 are constants. The term Be~^/^+ comes from differentiation of e~*'^^. 
This allows to solve for jS and 7 as a function of B. The renormalized cosmo logical constant 
A is defined as the variation of 7. From ( [4.5| ) we infer 

A = 7 - 7„it ~ Se-^/"+ (4.6) 

Next, we consider any correlation function ^-^^trA"^) (obtained by evaluating p{k)). It is 
clear from Eq. (|4.4| ) that its variation is of the form 

p{k) = C4 + C5Se-^/'^+ + C6e-^/'^+ + ■ ■ • (4.7) 

and therefore its singular part is 

P(^)sing ~ (4.8) 

This last result displays a zero string susceptibility exponent and a logarithmic correction 
characteristic of a c = 1 theory ||18|| . 



5. Summary of results and prospects 

We have given in this letter the exact solution of the six vertex model on random 
planar graphs at zero cosmological constant and described the vicinity of this critical line. 
Note that all the expressions found agree with the known results in the A = and A = 1 
cases. In particular, for A = the critical value 6crit = and the behaviors ( |3.1tj| ) 
coincide with earlier calculations , while for A = 1 the critical value 6crit = 4^ and some 
correlation functions such as (^j^trXX^'^ = ^(4 — tt) coincide with The limit A — 2 
is singular, as can be seen in ( |3.8| ), and therefore we cannot compare directly our results 



with those of ||T^,0; in fact, even though the A = 2 probably does have a critical point 



for a finite value 6crit [01 ' believed to belong to the universality class of pure gravity. 



which is different from the critical behavior found for A < 2 (c = 1 theories). 

The solution is explicit enough to allow us to give exact expressions for correlation 
functions at criticality, and the method can be used to generate an exact expansion around 
the critical region. Two critical exponents have been computed this way: the exponent 
governing the singularity of the density of eigenvalues, and the string susceptibility expo- 
nent; the latter turned out to be zero plus a logarithmic scaling violation, confirming the 
central charge c = 1 of the infra-red CFT. 

This raises the hope that it is possible to solve much more general models. For 
example, one can replace the quadratic potential of A with a more general one. This could 
be used to simulate e.g. vortices in the 6- vertex model or dilution in the deformed 0{2) 
model. There is in principle no problem to finding the critical properties of these models, 
and we expect a rich structure of multi-critical points. All this should appear in a future 
publication, as well as a more detailed description of the off-critical region. 

Finally, let us mention that these results should have an interesting application to 



knot theory and more precisely the counting of alternating links, see [|T9[. This is currently 
under study. 
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